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A REVIEW OF FLIGHT AND WIND TUNNEL MEASUREMENTS 
OF BOUNDARY LAYER PRESSURE FLUCTUATIONS 
AND INDUCED STRUCTURAL RESPONSE 

by David Alan Bies 
Bolt Beranek and Newman Inc. 


SUMMARY 


A review is presented of available data on boundary 
layer pressure fluctuations and induced structural response, 
from flight and wind tunnel investigations. The wind tunnel 
data include flat plate pressure fluctuation spectra and 
space-time correlations, displacement and acceleration spectra 
of flat flexible panels, and sound power spectra radiated 
by flat flexible panels. The flight data include pressure 
fluctuation spectra and "equivalent acoustic spectra, " the 
acoustic fields that would produce the same response as 
the aerodynamic fields. 

In order to use the same normalization procedure with 
all the data, engineering curves have been derived for 
estimating boundary layer parameters. These curves extend 
the estimates to Mach numbers up to 4 and Reynolds numbers 

q 

based on a characteristic length up to 2 x 1CK. 

The pressure fluctuation data show considerable 
scatter, especially in the wind tunnel investigations. 

The experimental results suggest that the scatter may 
be due to highly localized flow perturbations. It is 
argued that these perturbations may not be significant 
in determining structural response. General recommendations 
are given to guide experimental studies, A simplified 
procedure for estimating boundary layer pressure spectra 
is given in an appendix. 


INTRODUCTION 


With the advent of high performance flight vehicles, 
aerodynamically-induced vibration has become of great 
significance to the vehicle designers. Aerodynamic 


excitation may come about through various causes, such 
as flow over cavities or around steps, base pressure 
fluctuations, or loosely attached shock waves. One form 
of significance which has received much attention is 
that due to pressure fluctuations in the turbulent boundary 
layer over the vehicle skin. This source may be of parti- 
cular significance to the supersonic transport designer. 

A large literature exists that reports measurements 
of boundary layer pressure fluctuations over flat plates 
in wind tunnels and on the surfaces of various flight 
vehicles, evidencing a lively interest in such phenomena. 
Unfortunately, it is difficult to draw general conclusions 
from this literature because there is no uniformity of 
reporting the data. The interest which has prompted 
the many investigations has also prompted the writing 
of this report in which an attempt is made to intercompare 
published data and draw general conclusions. 

In most cases of reported flight vehicle measurements, 
no information is given about the boundary layer. If one 
is to compare boundary layer pressure fluctuation measurements 
with each other or with wind tunnel measurements, some 
estimates of the boundary layer parameters must be provided. 

In addition many wind tunnel measurements are reported 
in such a way that intercomparison is rendered difficult, 
and here too some means for estimating missing information 
is required. 

The path followed in this report follows from the 
stated objective and the problems listed above. Thus 
the first task is to construct a set of procedures and 
engineering curves by which boundary layer parameters 
can be estimated. This information is contained in the 
section labeled Boundary Layer Parameters, The literature 
and the procedures involved in generating the data given 
in this section are discussed in Appendix A. 

A second p ro blem arises in t he consideration, of .. the 
measurement of boundary layer pressure fluctuations due 
to the natural loss in sensitivity of a finite-size pressure 
transducer at high frequencies and short wavelength. 

This problem has been the subject of much discussion and 
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has been dealt with in various ways by various investigators. 
An effort has been made to give this matter consistent 
treatment in this report through work outlined in the 
section Microphone Size Correction. 

In the section on Wind Tunnel Measurements, the 
results of pressure fluctuation measurements and correlation 
studies of the turbulent boundary layer over flat plates 
are discussed. Based on this material, a simplified procedure 
for estimating boundary layer pressure spectra is given in 
Appendix B. In addition a number of published papers on 
wind tunnel investigations of the response of panels to 
turbulent boundary layer excitation are reviewed, and the 
principal findings are summarized. 

The results of a number of flight vehicle measurements 
of turbulent boundary layer pressure fluctuations and 
of response of structures to such excitation are given 
in the section on Plight Measurements. It is shown in 
this section that the available data form a fairly consistent 
and comprehensive picture of the phenomena. 

The last section called Conclusions and Recommendations 
summarizes the principal findings of this report and makes 
some general recommendations for further work. 


SYMBOLS 


C^ spatial part of the correlation function 

C 2 time part of the correlation function 

Co,Ch spatial part of the correlation function in long! 
^ tudinal and traverse directions 

c f local coefficient of skin friction 

c f local coefficient of skin friction in transformed 

coordinates (see ref. A7) 

D microphone diameter 

F( f ) pressure spectral density 

M local Mach number 
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free stream Mach number 

CO 

q free stream dynamic pressure 

R space-time correlation function 

r recovery factor 

Re Reynolds number based on distance from the wall 
to the sublayer (see ref. A7) 

Re Reynolds number based on distance x 

Re Reynolds number based on distance from wall 

J 

Re s Reynolds number based on boundary layer thickness 

Re^ Reynolds number based on momentum thickness 

Re^ Reynolds number based on momentum thickness in 

transformed coordinates (see ref. A7) 

t time 

T temperature 

T temperature in the sublayer (see ref. A7) 

s 

T w temperature at the wall 

free stream temperature 
T qoq free stream stagnation temperature 

U local flow velocity 

U convection velocity defined by means of the turbulent 

c boundary layer correlation properties 

free stream flow velocity 

x coordinate in plane of wall in direction of flow, 

distance from inception of turbulence on a flat plate 

y coordinate normal to the wall into the stream 
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z coordinate in plane of wall transverse to direction 

of flow 

a an exponent 

y ratio of specific heats 

5 boundary layer thickness 
_* 

0 displacement thickness 

R displacement transverse to direction of flow 

constant displacement transverse to direction of 
flow 

6 momentum thickness 

M- local gas viscosity 

free stream gas viscosity 

1 displacement parallel to direction of flow 

constant displacement parallel to direction of flow 

p local gas density 

p m free stream gas density 

t time delay 

if a function (equation Al8) 

cd frequency in radians per second 


BOUNDARY LAYER PARAMETERS 


As mentioned in the introduction, it seldom happens 
that a complete description of the boundary layer is 
presented with measurements of the boundary layer pressure 
fluctuations. This shortcoming is especially true of 
flight measurements and is also quite often evident in 
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published wind tunnel measurements. To remedy the total 
or partial lack of such information, engineering curves 
have been developed and are given in figures 1 through 4. 
The information given in these curves is based on empirical 
data discussed in Appendix A and is valid strictly for 
unperturbed flow over a flat plate. We shall use the 
information contained in the figures to supply missing 
parameters in order to compare published and unpublished 
boundary layer pressure fluctuation data. In doing so, 
we shall attempt to represent the actual measurement 
configuration by an equivalent flat plate model. 

The procedure which has been used in the preparation 
of this report is as follows. Typically, a report of 
flight measurements of boundary layer pressure fluctuations 
will give information about the location of the transducer, 
air speed, and altitude. Prom the information about the 
transducer location, an estimate of the characteristic 
length to the point of inception of turbulent flow on 
an equivalent flat plate is made. This characteristic 
length is quite often taken as the distance to some leading 
edge. With the altitude information and the characteristic 
length estimate just determined, an effective Reynolds 
number is computed from ARDC standard atmosphere tables 
(ref, l). From the given air speed and the altitude 
information, a Mach number is computed, again by use of 
the ARDC standard atmosphere tables. The three bits of 
derived information -■* characteristic length, Reynolds 
number, and Mach number -*■ are sufficient to determine 
the four boundary layer parameters — boundary layer 
thickness using figure 1, displacement thickness using 
figure 2, momentum thickness using figure 3* and local 
skin friction coefficient using figure 4. Our normaliza- 
tion procedure makes use of the displacement thickness 
and the local skin friction coefficient. 


MICROPHONE SIZE CORRECTION 


Experimental and theoretical Investigations have 
shown that the finite size of a pressure transducer limits 
its response in the upper frequency range. This happens 


6 


because frequency and spatial extent are related so that 
high frequency corresponds to small spatial extent. 

Thus, at some frequency determined by the ratio of micro- 
phone diameter to boundary layer displacement thickness 

D/5, cancellation between positive and negative pressures 
begins to occur over the transducer surface. 

This problem has been investigated by Corcos (ref. 2), 
Willmarth and Roos (ref. 3) and Chandiramani (ref. 4). 
Willmarth and Roos have shown theoretically and experi- 
mentally that the earlier work of Corcos overestimates 

the size correction for values of D/5* less than unity, 
but for values of unity or greater Corcos 1 s estimates 
are correct. Corcos* s estimates have been found to be 
in satisfactory agreement with the results of experiment 

for a value of D/6 much larger than unity, according 
to Gilchrist and Strawderman (ref. 3). Unfortunately, 
the latter authors have not specified the size of this 
parameter, but it may be estimated to be of the order 
of eight by reference to the work of one of their co-workers 
(ref. 6). 


Chandiramani has independently estimated the micro- 
phone size effect and has obtained slightly different 
results from those of Willmarth and Roos. The expected 
reduction in pressure sensitivity of a transducer due 
to its finite size according to these authors is shown 

in figure 3. The curve for D/6 of the order of unity 
or greater is due to Corcos. The family of solid lines 

for D/6 of 0.1, 0.2 and 0,5 is due to Wlllmarth and 

Roos. The family of two dashed lines for D/6* of 0.33 
and 0.66 is due to Chandiramani 

In the presentation of data Included In this report, 
microphone corrections according to Corcos, Willmarth 
and Roos have been Included where sufficient Information 
was given by the author. However, in some cases It is 
not clear that some such corrections have not already 
been made. In cases of doubt no corrections have been 
made. If the data has been corrected, this fact Is indicated 
In the corresponding figure. 
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WIND TUNNEL MEASUREMENTS 


Boundary Layer Pressure Fluctuations 

Boundary layer pressure fluctuations have been measured 
fey many investigators. Those whose work has feeen included 
in this review (refs. 7 through 19) are listed in Table I. 

For convenience the references have feeen listed roughly 
in order of increasing Mach number, and some of the parameters 
pertinent to particular Investigations are also listed. 

The table serves as an index to the wind tunnel measurements 
presented In figures 6 through 19. 

The investigators listed in Table I have reported 
their results in many different ways, so that it 13 difficult 
to determine a general concensus. Following procedures 
described earlier, we have attempted to determine reasonable 
estimates for quantities which are sometimes missing in 
published reports (such as the local coefficient of skin 
friction), and we have used such estimates to recast 
published work in the form shown in figures 6 through 19. 

The normalization presented in these figures or very similar 
normalizations have been used fey a number of investigators, 
and In those cases little or no change has been made in 
the data as originally published. 

The first thing that comes to mind as one scans 
the data shown in the figures is the question. How do 
the various measurements compare?" Such a comparison 
is shown in figure 20. Below a Strouhal number of unity, 
there Is quite a bit of agreement among most experimenters 
on an upper bound of the spectrum. Above a Strouhal 
number of unity, the upper bound of the spread is determined 
by the work of Kis tier and Chen (ref. 17). The lower 
bound for all Strouhal numbers is determined by the work of 
Harrison (ref. 8). In the latter case it Is not clear 
whether or not corrections for finite microphone size 
were included by Harrison, so that no corrections have 
been made In this report. 

As an aid to the reader a line has been placed where 
most of the data seem to cluster. Where the line Is 
discontinued, the data Is about equally distributed over 
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TABLE I 


MEASUREMENTS OF BOUNDARY LAYER PRESSURE FLUCTUATIONS 


Investigator 

Tunnel 
Size + Type 

Mach 

No. 

Reynolds 

No. 


Fig. 

No. 

Jorgensen. Donald W. 
(ref. 7} 

Continuous 
18 " x 15" 

0.04 

0.18 

(x)5.9 x 105 
2.4 x 10° 

6.7 

11.4 

6 

Harrison. M. 
{ref. 8) 

Continuous 
15" x 20" 

0.04 

0.18 

(x ) 1 X lOg 

4 x 10° 

1.2 

7 

Baroudi, Ludwig, 

+ Ribner (ref. 9) 

Continuous 
8" x 12" 

0.04 

0.15 


0.092 

8 

Bakewell, H. P. Jr., 
et al (ref. 6) 

Continuous 
3-1/2" diam. 

0.1 

(x)4 x 10° 

0.31 

9 

Wilmarth + Wooldridge 
(ref. 10) 

Continuous 
60 " x 84" 

0.2 

(x)3.1 x lol 
3.8 x 10' 

4.0 

10 

Bull, Wilby + 

Blackman (ref. 11) 

Continuous 
9" x 6" 

0.3 

0.5 

(x)0,5 x loZ 

3 x 10' 

0.15 

0.51 

11 

Murphy, J. S., et al 
(ref. 12) 

Blowdown 
12" x 12" 

0.59 

(0)1.3 x 10 4 

2.77 

12 

Serafini, J. S. 
(ref. 13) 

Continuous 
8" x 18 " 

0.6 

(x)l X lOo 
1.4 x 10 8 

0.195 

1.33 

13 

Maestrello, L. 

(refs. l4 and 15 ) 

Continuous 
7.3" x 3.5" 

0.63 

0.78 


0.45 

14 

Williams, D. J. M. 
(ref. 16) 

Continuous 
6" x 2-1/2" 

1.2 

1.5 


5 

8 

15 

Kistler + Chen 
(ref. 17) 

Continuous 
18 " x 20" 

3.5 

(9)1.55 x 10 q 

1.1 

16 

Murphy, J. S., et al 
(ref. 12) 

Blowdown 
12" x 12" 

3.46 

(0)4,9 x 10^ 

0.43 

17 

Westley, R. 
(ref. 18 ) 

Blowdown 
5" x 5" 

2.91 

3.96 
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Speaker + Allman 
(ref. 19) 

Blowdown 
12" x 12" 

0.42 

3.45 

(9)1.6 x 107 
5.3 x 10 

0.17 

1.3 

19 


(x) Reynolds number based on distance from inception of turbulence. 
(9) Reynolds number based on momentum thickness. 
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the indicated range. The large spread leads one to ask 
if a better normalizing scheme exists. For the consideration 
of this question the work of Speaker and Ailman is useful 
(ref. 19). 

Speaker and Ailman report a large number of measurements 
of unperturbed flow over a flat plate for a large range 
of Mach numbers. In addition they provide a rather complete 
list of various boundary layer parameters. Their report 
shows quite a large variation in measured pressure fluctua- 
tions across the plate, for any one Mach number. It is 
clear that no choice of the given parameters would bring 
their data together for a single Mach number, since none 
of the parameters vary significantly across the plate. 

However, Speaker and Ailman do present a figure which 
summarizes the mean of their findings for each Mach number. 
These mean values are shown in figure 19. 

Figure 19 shows the following interesting trend. 

Above a Strouhal number of two where all of the data 
curves roll off sharply, the trend is for the low and 
high Mach number measurements to lie furthest to the right, 
and the measurements close to Mach number unity furthest 
to the left. A review of the boundary layer parameters 
presented by Speaker and Ailman shows that the boundary 
layer thickness, displacement thickness, and momentum 
thickness show the same general behavior, being largest 
at lowest and highest Mach numbers and smallest near unity 
Mach number. Thus, normalization using either of the other 
boundary layer thicknesses in computing the Strouhal number 
would show the same general trend as shown by normalization 
based on displacement thickness. 

One might wonder if some combination of parameters, 
perhaps Involving a Reynolds number based on one of the 
thicknesses, might give better normalization. The author 
has not found any such dimensionless combination of parameters 
which could significantly coalesce the high-frequency 
rolloff portions of the spectra. It is clear that the 
behavior shown by the data of Speaker and Allman follows 
from the choice of normalization parameter, displacement 
thickness, but it is not clear what other choice or combina- 
tion of parameters should be made. Figure 20 unfortunately 
appears to represent the present state of the art for 
wind tunnel measurements. 


Space-Time Pressure Correlations 

Interest in the space-time surface-pressure correlation 
of the turbulent boundary layer stems primarily from the 
fact that a knowledge of the complete correlation is suffi- 
cient to calculate the response of a structure excited 
by the boundary layer. It has been pointed out (ref. 4) 
that the response of a structure excited by a turbulent 
boundary layer is unfortunately rather sensitive to the 
fine details of the space-time surface-pressure correlation. 

In addition to the primary interest In pressure correlations 
for Investigating structural response, correlation measurements 
are of interest to provide Insight into the basic mechanics 
of the turbulent boundary layer. For example, estimates 
of the convection velocities (refs. 10, 11, 13, 19, 20, 21) 
lifetimes (refs. 10, 22) of the eddies which produce the 
random pressure fluctuations in the turbulent boundary 
layer have been deduced from correlation measurements. 

Table II presents a summary of wind tunnel pressure 
correlation measurements. A review of these measurements 
indicates that, although the investigators have obtained 
quite detailed measurements of the correlation properties, 
there exists disagreement In the different measurements 
of some of the basic properties of the surface-pressure 
correlations. Efforts to relate these experimental discrep- 
ancies to differences In the flow conditions have been 
frustrated by the present lack of knowledge of the basic 
physics of the turbulent boundary layer. 

Terminology . - Before discussing the wind tunnel 
measurements, it is convenient to define a coordinate 
system and review some terminology relevant to cross- 
correlations. Let us define an orthogonal coordinate 
system x, z In the plane of the wind tunnel wall. Let 
x be the coordinate in the direction of flow and z be 
the coordinate In the direction transverse to the flow. 

The complete space-time surface-pressure correlation 
R(x, z, t; x + 4, z + T), t + t) is a croes-correlation 
obtained by multiplying together two surface-pressure 
measurements and then time -averaging the resultant product. 
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TABLE II 


SPACE-TIME CORRELATION MEASUREMENTS 
OF BOUNDARY LAYER PRESSURE FLUCTUATIONS 


INVESTIGATOR 

Bull, Wilby 
and 

Blackman 
(ref. 11 ) 

Speaker 

and 

Allman 
(ref. 19) 

Serafinl 
(ref. 13 ) 

MACH NUMBERS 

' ' 

0.3 and 0.5 



0.59, 1.41, 
3.45 and 3.48 

0.6 

5* (in.) 

approx . 

approx. 

0.32 

0.05 to 0.15 

0.03 to 0.35 


, * 
D /6 

approx. 

approx . 

0.38 


0.15 to 0.51 

0.16 to 2.6 


MEASUREMENTS REPORTED 




1 . Autocorrelations 

X 


X 

2. Longitudinal 

X 

X 

X 

correlations 




3. Transverse 

X 

X 

X 

correlations 




4. Off-Axis 

X 

X 

X 

correlations 




5 . Spatial correlation 

X 

X 


contours 




6 . Convection velocity 

X 

X 

X 

vs longitudinal 
sepration 




7. Convection velocity 

X 



vs frequency 





Correlations 

Correlations 

Correlations 


of narrow- 

in perturbed 

of low and 

SPECIAL TYPES OF 

band pressure 

supersonic 

high pass 

MEASUREMENTS 

fluctuation 

flow 

filtered 

measurements 

Effects of 
assuming 
separable 
correlations 

pressure 

fluctuation 

measurements 
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TABLE II (cont.) 


SPACE- TIME CORRELATION MEASUREMENTS 
OF BOUNDARY LAYER PRESSURE FLUCTUATIONS 


INVESTIGATOR 

Maestrello 
(refs. 14, 15, 22 , 23 ) 

Willmarth 

and 

Wooldridge 
(ref. 10 ) 

MACH NUMBER 

0.52 

approx . 
0.18 

5* (in.) 

0.155 

0.5 

D/5 

0.4 

0.33 

MEASUREMENTS REPORTED 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

1. Autocorrelations 

2 . Longitudinal 
correlations 

3. Transverse 
correlations 

4. Off-axis 
correlations 

5. Spatial correlation 
contours 

6 . Convection velocity 
vs longitudinal 
separation 

7. Convection velocity 
vs frequency 

SPECIAL TYPES OF 
MEASUREMENTS 

Correlations 
of narrow- 
band pressure 
fluctuation 
measurements 
Eddy lifetime 
vs Mach number 

Spatial 
decay of 
wave number 
components 
3-D plot of 
longitudinal 
space-time 
correlation 
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One pressure measurement is taken at the point (x, z ) 
on the surface at time t , and the other measurement is 
taken at a different point (x + i, z + r) on the surface 
at a different time T + t. 

In all the investigations summarized in Table II, 
except the perturbed flow study of Speaker an ^ 

(ref. 19), it is assumed the surface-pressure field is 
spatially and temporally homogeneous, so that the correlation 
is a function only of the spatial differences R, and 
the temporal difference r 9 and can be written as R(s, R* 

The inhomogeneous cases considered in reference 19 give 
rise to correlations which are functions of the measurement 
positions as well as the separations. In actuality the 
correlations reported by the various investigators are 
special cases of the complete space-time surface-pressure 
correlation. We have used the following terminology 
in Table II and In the text: R(0, 0, t) autocorrelation; 

r(£ 0, 0) longitudinal correlation; R(0, R, 0) transverse 

correlation; contours of constant R(l, r, 0) “ spatial 
correlation contours; R 0, t) longitudinal space-time 

correlation (£ Q indicates a constant spatial separation); 

and R(0, r 0> t ) transverse space-time correlation. 

Correlation results . - The cross correlation is 
the Fourier transform of the cross power spectral density. 

The autocorrelation Is the transform of the power spectral 
density function. This relationship has been used by 
Speaker and Allman (ref. 19) to determine the power spectral 
density of the turbulent boundary layer pressure fluctuations. 
The same relationship has been used by others (refs. 11 
and 22) as a check on their correlation measuring equipment. 

Measurements of the longitudinal and transverse corre- 
lations reported by the authors listed in Table II are » 

presented in figures21 and 22. With one exception (Maestrell ) 
the subsonic data in figure 21 show that the correlation 
function drops to zero in the direction of flow within 
a space of the order of three to five times the displacement 
thickness, and stays negative. The correlation gradual y 
approaches zero at some point beyond twenty times the 
displacement thickness of the boundary layer. T ^ e cme 
supersonic measurement Indicates the same general trend 
but the zero crossing occurs in a much shorter distanc . 
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In plotting the data of Willmarth and Wooldridge 
(ref. 10} in figure 21 note has been taken of the fact 
that their space-time correlations show the presence of 
an acoustic disturbance traveling upstream. This disturbance 
had the effect of introducing positive correlation near 
zero time delay which only became significant for values 
of the transducer separation distance of the order of 
5 to 7 times the displacement thickness. In replotting 
the data the hump in the correlation function due to the 
acoustic signal has been edited out in order to present 
the estimate of the longitudinal correlation shown in 
the figure. 

The data of figure 12 show that the correlation 
laterally remains positive and gradually approaches zero 
at a distance greater than twenty times the boundary layer 
displacement thickness. Again the one piece of supersonic 
data indicates that the correlation approaches zero in 
a shorter distance than does the subsonic data. The 
trends indicated in figures 21 and 22 have been observed 
by other workers as well (for example, see refs. 6 and 24). 

As shown in Table II, In addition to longitudinal 
and lateral spatial correlation, Investigations have been 
carried out at various off-axis directions to the flow. 
Several authors have also presented contour plots of 
constant correlation R(£, tj, 0) (refs. 11, 19, 21). 

As might be expected from Inspection of figures 21 and 
22, such plots strongly show the differences between the 
various experimental results and do not look at all alike. 

It is suggested that one might draw mean curves through 
the data given in figures 21 and 22 and use these to obtain 
estimates of contours of constant correlation. The off-axis 
spatial correlations indicate that such contours would take 
on a figure-eight configuration. 

References 10, 11, 13, 19, and 22 listed in Table II 
all contain measurements of the longitudinal space-time 
correlation R(s , 0, t) 9 with a fixed longitudinal transducer 
separation ^References 11, 13, and 19 also contain 

measurements of the transverse space-time correlation 
R (0, %, T ) obtained for fixed transverse transducer 
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separation. In addition, reference 11 contains longitudinal 
and transverse space-time correlations of surface-pressure 
measurements which have been filtered in narrow bands 
and reference 10 contains longitudinal and transverse 
space-time correlations of surface-pressure measurements 
which have been filtered by high-pass and by low-pass 
filters. 

Speaker and Ailman (ref, 19) considered the possibility 
of factoring complete space- time correlation into three 
functions, each of which depends only °n one sep a ration 
variable. For the flow conditions investigated, it was 
concluded that the correlation could be expressed as the 
product \ (I, r») C 2 (r) of one function depending only 

on the spatial separation variables and another function 
depending only on the time delay. However, it was concluded 
that separation of the spatial correlation function C 1 (4,tj) 

into two functions C^{i) and C^(r|) resulted in considerable 

error even for the case of unperturbed supersonic flow, 
in addition, it was concluded that, in the case of supersonic 
flows nerturbed by steps or shocks, the spatial P al *p 
c(!!nfoftheco?r elation behaved ’differently in different 

«m»tlal r^eions of the flows, whereas the temporal part 
C P tr) of ?he Correlation behived similarly in the different 

spatial regions of the investigated flows. 

Convection Velocities . - Figure 23 presents estimates 
of the broad band convection velocity as a function of 
longitudinal separation determined by various authors 
listed in Table II. [The convection velocity corresponding 
to the longitudinal separation is usually calculated 

from measurements of the longitudinal space-time correlation 
R(| , 0, r) by dividing the separation by the time 

delay t for which the correlation R(L, °> T ) a maximum. 

Tt has been argued (ref, 4) that a more acceptable, b 
more difficult, method of measuring the convection velocity 
would be to divide the separation 1 at which the l°£S dtudinal 
space-tlme correlation R(l, 0, t q ) Is a maximum by th 

fixed time delay t q ,] The convection velocities shown 
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in the figure all start at approximately 0.5 (one 

half the free stream velocity), increase quite rapidly 
with increasing longitudinal separation, and level off 

at a value of approximately 0.8 U at a longitudinal 

°° * 

separation of approximately 10-30 6 . Since experimental 
measurements (refs. 10, 22) indicate that the high frequency 
eddies decay more rapidly than the low frequency eddies, 
the increase in convection velocity with increasing longi- 
tudinal separation is usually interpreted as evidence 
that the high frequency eddies travel slower than low 
frequency eddies. 

Bull Wilby and Blackman (ref. 11), Bakewell (ref. 6) 
and Corcos (ref. 21) have investigated the convection 
velocity as a function of frequency (dispersion) obtained 
from narrow-band measurements of the longitudinal space - 
time correlation. These measurements are presented in 
figure 24, where the dimensionless convection velocity 
has been plotted as a function of dimensionless wavelength 
(reciprocal frequency). The measurements of Bull et al 
have been taken directly from the reference, but in order 
to include the work of Bakewell and Corcos, an estimate 
of the displacement thickness was required. In both 
cases fully turbulent flow in a pipe was investigated. 

We assumed that the displacement thickness was 1/1 6 of the 
pipe diameter. 

The measurements shown in figure 24 indicate that 
the high frequency, short wavelength components travel 
more slowly than the low frequency long wavelength components. 
The measurements are consistent with the physical picture 
that the high frequency components are produced by small 
scale eddies, which are located close to the surface 
and therefore do not travel as fast as the low-frequency 
large-scale eddies, which project further out into the 
boundary layer. 

The data in figure 24 may reasonably well be represented 
by a straight line with a slope of the order of 0.09. 
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Structural Response 


Work on structural response to turbulent boundary- 
layer excitation has been carried on by Lucio Maestrello 
of the Boeing Company (refs. 14, 15, 22, 23), J. p * 
et al of the University of Southampton (refs. 11, 25, 2b ) 

G. R. Ludwig of the University of Toronto (ref. 24). 

D. A. Bies of Bolt Beranek and Newman Inc. (ref. 27), 
and by Paul F, R. Weyers of California Institute of Technology 
(ref. 28). The first four investigators placed a test 
panel flush with the side wall of a wind tunnel for their 
investigation. The last investigator listed above investi- 
gated the response of a thin Mylar tube to turbulent flow 
through it. Information pertinent to these tests is given 
in Table III. 

In all cases of response Investigation except that 
of Bies, the displacement response was measured, and the 
resulting spectra were dominated by the fundamental modes 
of the panel or tube and the first few overtones. Workers 
at Southampton have identified the various panel modes 
and have compared measured frequency and amplitude with 
the predictions of theory outlined by them. The theory 
predicts fairly well the resonant frequencies observed 
but gives poor agreement with the measured amplitudes. 


The workers at Southampton have used their theory 
to predict trends in response with variation of various 
parameters such as boundary layer displacement thickness 
and panel thickness. They have then sought to verify 
experimentally the dependence of panel response on these 
parameters, with mixed success. 


In addition to Investigating the response of panels 
to boundary layer excitation, the workers at Southampton 
have also investigated the response of their panels to 
a siren and an air jet. Their results indicate that the 
acoustic noise sources and the turbulent boundary layer 
are equally effective in producing displacement for a 
unit mean square pressure fluctuation. This matter is 
discussed further later In this section. 
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TABLE III 


WIND TUNNEL STRUCTURAL RESPONSE MEASUREMENTS 


INVESTIGATOR 

PANEL DIMENSIONS 

MODAL DENSITY 
RANGE 

MACH 

NUMBER 

DISPLACEMENT 

THICKNESS 

(inch) 

QUANTITY j 

MEASURED 

Size 

(Inch) 

Thickness 

(inch) 

Lucio Maestrello 
(refs. 14,15,22, 
23) 

aluminum 
7 x 12 

. 

I 

fundamental 
and first few 
overtones 

0.34 

to 

0.64 

0.155 

radiated power, 
displacement 

J. F. Wilby 
Bull, Wilby and 
Blackman 

Wilby and Richards 
(refs. 4,25,26) 

steel 
4x2 
3.5 x 3.5 
4 x 2.75 

0.005 

0.010 

0.015 

fundamental 
and first few 
overtones 

0.3 

0.5 

0.053 

to 

0.176 

displacement 

G. R. Ludwig 
(ref. 24) 

steel 
11 x 11 

0.0015 

0.002 

0.004 

0.008 

fundamental 
and first few 
overtones 

0.05 

to 

0.17 

0.054 

0.43 

radiated power, 
displacement 

Paul F. R. Weyers 
(ref. 28) 

mylar 
cylinder 
1.0 dia. 

0.005 

0.001 

0.0021 

fundamental 
and first few 
overtones 

0.11 

to 

0.28 


radiated power 

D. A. Bies 
(ref. 27) 

steel 
11.5 dia. 

0.012 

0.022 

high modal 
density 


0.068 

0.306 

acceleration 

































Maestrello and Ludwig have Investigated both panel 
response and radiated sound power from panels of comparable 
size. As Table III shows, however, Ludwig has looked at 
very thin panels, whereas Maestrello's panels are an 
order of magnitude thicker. In both of these cases the 
panel response was dominated by the fundamental mode of 
the panel and its first few overtones, Ludwig concentrated 
his effort on the measurement of radiated sound power 
and reports only one measurement of panel displacement. 
Maestrello reports many measurements of both displacement 
and radiated sound power. 

Typical of the spectra of the radiated sound power 
observed by Maestrello and also of Ludwig and Weyers is 
the set of curves shown in figure 25. In the figure is 
shown the sound power in one- third octave bands of noise 
radiated by a 0.080 Inch thick panel Into a reverberant 
room. The peaks in the spectra may be associated with 
the fundamental frequency of the panel and its first few 
overtones , 

From the figure we notice that two things happen 
as the flow velocity increases: The general level rises, 

and the response extends to higher frequencies. The first 
observation is understandable in terms of Increased pressure 
fluctuation levels with Increased flow, since the form 
of the pressure normalization parameter Implies that the 
low-frequency pressure spectrum increases proportionally 
to the cube of the flow velocity. The second observation 
may be understood on the following basis. If we define 
a hydrodynamic coincidence frequency as the frequency 
at which a flexural wave in the panel is equal to the 
mean convection speed of the turbulent boundary layer 
(taken as about 0.8 of the free stream velocity), we find 
that this frequency is directly proportional to the square 
of the flow velocity and Inversely proportional to the 
panel thickness. The spectra shown in the figure begin 
their sharp upper frequency roll off near this coincidence 
frequency in each case. Thus, the strongly excited portion 
of the panel response spectrum lies between the fundamental 
frequency and the coincidence frequency. 
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Maestrello, Ludwig, and Weyers all report that the 
radiated sound power is proportional to the fifth power 
of the free stream velocity. This power law dependence 
has been observed in at least one flight measurement 
(ref. 29). This relationship is not precisely followed, 
however, since power law variations about the fifth power 
were also observed. In fact, Maestrello showed that, at 
high Mach numbers for all but his thickest panel, the 
power relationship abruptly decreased to the 2.5 power 
law. Data which illustrates these findings taken from 
Maestrello’s work is shown in figure 26. 

The dependence of radiated power or panel displacement 
response on panel thickness may be inferred from the data 
in figure 26. Maestrello, Ludwig, and Wilby found that 
the panel displacement is proportional to some inverse 
power of the panel thickness. The power ranges between 
1.0 and 1.6. Maestrello summarizes these findings and 
compares them with the theoretical predictions of various 
workers whose predictions range between 1.0 and 3.0. 

Figure 26 illustrates another fact discovered by 
these investigations. The radiated sound power follows 
the mean square amplitude of response of the panels when 
the low order modes dominate the response. This conclusion 
is borne out also by comparisons of displacement and 
radiated power spectra given by Ludwig and Maestrello. 

This conclusion is contrary to what would be expected if 
the radiation impedance of the panel had a significant large 
real part, since in this case the radiated power should 
follow the panel velocity rather than the displacement. 

It should be remembered that if the velocity remains 
constant with increasing frequency, the displacement should 
decrease in proportion to the frequency, so that the effect 
discussed here should be readily observable. 

An effort was made to find relationships which would 
allow one to scale Maestrello’s and Ludwig* s work, but 
no simple scaling laws could be found. Ludwig gives 
relationships which seem to hold for his data, but they 
cannot be used to scale. Maestrello* s results. This diffi- 
culty may come about because Ludwig’s panels may be controlled 
by membrane forces rather than panel stiffness forces. 

This point is not clear. 
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In summary we may say that the investigations of the 
first four investigators listed in Table III have shown 
the following: 

(1) The displacement response of a lightly damped 
panel to turbulent boundary layer excitation 
is always dominated by the fundamental of the 
panel and its first few overtones. 

(2) The frequency range most strongly excited extends 
up to the frequency of coincidence between flexural 
waves and the convected turbulent boundary layer. 
Above this frequency the displacement spectra 

fall off rapidly. 

(3) In the frequency range of low modal density 
the radiated sound power is proportional to 
the mean panel displacement rather than panel 
velocity. 

(i|) In the frequency range of low modal density 
the radiated sound power Is proportional to 
the fifth power of the free stream velocity 
over a range of low Mach numbers, but at high 
subsonic Mach numbers the power law decreases 
to about the 2.5 power. 

(5) No theory apparently exists which gives good 
quantitative agreement between predicted and 
observed response, and no empirical relationships 
have been discovered by which one set of data 
may be scaled to that of another set of data. 

As mentioned earlier and as pointed out in Table III, 
Bies has carried on quite a different investigation of 
panel response. In his experiment the panels investigated 
were highly damped, and the acceleration response (rather 
than displacement response) was measured. Consequently, 
the low order modes were deliberately suppressed. The 
frequency range of high modal density was investigated, 
where it was shown that fairly good agreement between 
theory and experiment exists. In actual flight measurements 
the dominant portions of the internally observed spectra 
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(for example, the vibration spectra from which curves 
of equivalent acoustic fields have been generated as 
in figures 36 through 4l) are in the frequency range 
of high modal density for the skin panels. Thus, Bies 
claims that his results may be scaled to estimate response 
of full-scale supersonic transport panels. 

In addition to investigating response of panels 
to unperturbed flow, Bies has also Investigated the effect 
of mild perturbations such as an Impinging shock and an 
expansion. He has also investigated the response of his 
panels to a reverberant acoustic field and has compared 
the response to such excitation to the response to a 
turbulent boundary layer. In addition, Bies's Investigation 
has been carried out in the supersonic flow range, while 
all other investigations have been In the subsonic flow 
range . 

The following conclusions are drawn from Bies’s report: 

(1) Fairly good quantitative agreement was obtained 
between theoretically predicted and measured 
panel response to turbulent boundary layer 
pressure fluctuations for the cases of unperturbed 
flow. 

(2) To a fair approximation, one could associate 

a single transfer function with a panel, relating 
acceleration response to turbulent boundary 
layer pressure fluctuations under conditions 
of unperturbed or mildly perturbed flow. 

( 3 ) At Mach 3.5 boundary layer pressure fluctuations 
were of the order of 6 dB more effective in 
driving the 0 . 22 " thick panel than comparable 
reverberant acoustic fields. This Is In agreement 
with theory, which predicts the ratio to be 

7 to 8 dB. 

(4) At Mach 1,4 two regimes were observed. In the 
low-frequency range, a reverberant acoustic 
field was slightly more effective in driving 
the panels than were boundary layer pressure 
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fluctuations. In the high-frequency range, 
the boundary layer pressure fluctuations were 
as much as 7 or 8 dB more effective in driving 
the panel than the reverberant acoustic field. 
This observation is in qualitative agreement 
with the predictions of theory. 


This last result appears also to be in qualitative 
agreement with the Southampton work which showed equivalence 
between acoustic and aerodynamic excitations. The Southampton 
investigation was in the low-frequency low-Mach number 
range, where Bies’s work would indicate that such equivalence 

would exist. 


FLIGHT MEASUREMENTS 


Boundary Layer Pressure Fluctuation Measurements 

Flight measurements of boundary layer pressure fluctua- 
tions have been assembled from various sources (refs. 2y 
through kl) and are presented in figures 27 through 42. 

In most cases the data come from published reports, while 
In the cases of Titan II, Titan III, and Saturn the data 
have been acquired by the author through contact with the 
people who made the measurements. Microphone size correc- 
tions have been made by the author wherever it was clear 
that the corrections were necessary. 

In practically all cases of flight measurements 
considered here, no information about the boundary layer 
was available. In order to normalize the available data, 
the procedure outlined above was used. The altitude, 
air speed or Mach number, and a characteristic dimension 
were used to estimate a Reynolds number. From the < estimated 
Reynolds number, characteristic dimension, and Mach number, 
estimates were made of displacement thickness and coeffi- 
cient of skin friction. The normalizations used In most 
of the figures are thus necessarily approximate. However, 
the data contained in figures 27 and 30 are based on 
measured values of the boundary layer parameters. 
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The data in the figures have been placed in two 
groups, with measurements on airplanes in the first group 
and measurements on space vehicles in the second group. 

In general an effort was made to place the airplane 
microphones in regions of relatively smooth flow, but 
this was not possible for the space vehicle microphones. 

In fact, the space vehicle microphones were often in 
regions of highly perturbed flow (shocks or separated flow). 
Thus our figure grouping places the less perturbed flows 
first and the more perturbed flows last. Following this 
general line, the various measurements have been ordered 
in order of Increasing Mach number within each group, 
assuming that flow will more likely be perturbed at higher 
Mach numbers than at lower Mach numbers. 

The airplane flight data contained In figures 27 
through 34 show the following trends: a general buildup 

in the low-frequency portion of the normalized spectrum 
with Increase in Mach number, and a general asymptotic 
approach to the same values at high frequencies for all 
Mach numbers. The glider data contained in figure 27 
form a lower bound for low Mach number and least perturbed 
flow. The glider data also define the asymptotic behavior 
of all data at high frequencies. The C-141 data of figure 30 
form a kind of upper bound for relatively unperturbed flow 
but high subsonic Mach number at low frequencies. The 
data of figure 33 appear to be rather erratic for reasons 
unknown and perhaps should not be compared with the other 
data, since they were taken on a nose cone which does 
not simulate a flat plate very well. In summary, the 
airplane flight data appear to be fairly self-consistent 
and in line with wind tunnel measurements. 

The spectrum levels observed on various space vehicles 
shown in figures 35 through 4l are generally quite high 
in the low portion of the spectrum but are In line with 
the Idea that they might asymptotically approach the 
high-frequency spectrum levels of unperturbed flow of 
figure 27. The very high levels observed on the space 
vehicles are thought to be due to highly perturbed flow. 

The data shown In figure 42 were taken on a Mercury 
test vehicle. The configuration was such that the flow 
was quite strongly perturbed. The very high levels In 
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the low-frequency range are thought to be due to the 
presence of a spoiler which strongly affected the subsonic 
and just sonic flow but was relatively ineffective at 
supersonic Mach numbers. 

Structural Response 

In figures 32 and 3 6 through 4l an additional bit 
of information has been included. In these figures we 
have inserted solid lines labeled "equivalent acoustic 
field.’' The lines are based on measurements of strain 
or acceleration on the skin or of acceleration on skin- 
mounted components under two conditions. The structure 
is excited by acoustic noise in the first condition and 
by turbulent boundary layer pressure fluctuations in the 
second. The lines represent the acoustic spectrum which 
would produce the same structural response as observed 
when excited by the boundary layer pressure fluctuations. 
Comparison of the equivalent acoustic pressures (represented 
by the line) with the measured boundary layer pressure 
fluctuations (shown as points) shows the relative efficiency 
of acoustic and aerodynamic exicitation. 

For most of the measurements, the equivalent acoustic 
pressure spectra are close to the boundary layer spectra. 

From this we conclude that acoustic and aerodynamic pressure 
are approximately equally efficient in exciting structural 
response, for the configurations measured. However, the 
data shown in figures 37 and 38 are significant exceptions 
to this conclusion. In these two cases acoustic and 
aerodynamic exictation appear to be equally effective 
over the low-frequency portion of the spectrum, but above 
a Strouhal number of about 0.2, the acoustic excitation 
becomes less effective, as indicated by the higher equivalent 
acoustic levels required to achieve the same response as 
produced aerodynamically. This general behavior Is in 
line with Bies*s observations on panel response in a wind 
tunnel, discussed above. The order of magnitude difference 
in efficiency from the flight data is also In agreement 
with Bies*s results. 
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CONCLUSIONS AND RECOMMENDATIONS 


The principal conclusions of this report are as 
follows : 

(1) Fourteen published wind tunnel measurements 
of boundary layer pressure fluctuations show 
spectra that vary widely one from another. 

However, over a limited frequency range. 

* . 

0.05<u)5 /U co <5.0, sufficient experimental data 

agree to define a characteristic spectrum 
about which many of the measurements cluster. 

(2) Flight measurements of boundary layer pressure 
spectra in unperturbed flows are in general 
agreement with each other and with the charac- 
teristic spectrum of the wind tunnel measurements. 

( 3 ) Flight measurements of perturbed flow generally 
show a buildup of the low frequency portion 

of the spectrum. The spectrum asymptotically 
approaches the unperturbed flow spectrum at 
high frequencies. 

( 4 ) There is general agreement that the mean convection 
velocity determined by space-time correlation 
measurements of the turbulent boundary layer 
increases with spatial separation between obser- 
vation points and that the asymptotic value 

for large separation distances is about 0.8 
of the free stream velocity, 

( 5 ) The narrow band convection velocity is approximately 
proportional to (cub*/^)" for subsonic flow. 

(6) Longitudinal and lateral spatial correlation 
are in general agreement though they differ 
widely in detail. As a consequence there is 
a general lack of agreement between various 
investigators on the shape of the contours 

of constant correlation in the turbulent boundary 
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layer, although the measurements have all been 
made in approximately the same Mach number 
range. 

(7) Measurements of structural response of panels 
excited by boundary layer turbulence in wind 
tunnels have shown that the displacement response 
is dominated by the first few modes of the 
panel. At low flow velocities the displacement 
response is proportional approximately to the 
fifth power of the flow velocity, but at high 
flow velocities the dependence decreases to 

the 2.5 power. No general theory or empirical 
relations exists which enable the scaling of 
one set of panel displacement response measurements 
(first few modes) to another set of measurements. 

(8) Sound power radiated by aerodynamically excited 
panels is proportional to the mean square displace- 
ment of the panel and thus follows the power 

law relationships to flow velocity described 
in item 6. 

( 9 ) In the frequency range of high modal density, 
one set of acceleration response measurements 
shows reasonable agreement with the predictions 
of theory. At low frequencies acoustic and 
aerodynamic excitations are approximately equally 
effective in exciting panel response. At high 
frequencies aerodynamic excitation may be as 
much as 8 dB more effective than acoustic noise 
in exciting panel response. The frequency of 
crossover between these two ranges decreases 
with increasing Mach number, but the general 
conditions which determine this crossover frequency 
quantitatively remain to be determined, 

(10) Measurements of the structural response of flight 
vehicles show a relative response to acoustic 
and aerodynamic excitations in qualitative agree- 
ment with the wind tunnel measurements described 
in item 9. 
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(ll) Although there have been many measurements of 

surface pressure fluctuations, very few measure- 
ments include sufficient aerodynamic data to 
permit comparisons with other data, or to extrapo- 
late the data to other conditions. In order 
to generate the aerodynamic data, a system of 
engineering curves has been developed from 
published measurements of local skin friction 
and boundary layer profiles to enable the esti- 
mation of boundary layer parameters for Mach 
numbers up to 4 and Reynolds numbers based on 

Q 

a chacteristic length up to 2 x 10 . Use of 
this system of curves has enabled the uniform 
presentation of flight and wind tunnel surface 
pressure measurements for intercomparison In 
this report. 

In addition to the specific conclusions drawn from 
the body of material reviewed In this report, some general 
conclusions and recommendations may be drawn. A great 
many wind tunnel measurements of turbulent boundary layer 
pressure fluctuation have been reviewed. In spite of 
apparent care taken by the various investigators, their 
spectrum and correlation results show wide variation when 
compared with each other, although each investigator generally 
presents data that are quite self consistent. One exception 
exists, however, where the investigators made a great 
many measurements over an extended region in the test 
section of their tunnel (ref, 19). Their results showed 
rather wide scatter, but their results are well within 
the scatter of all of the measurements discussed in this 
report. One wonders whether other experimenters have 
looked over extended regions of their test sections and 
found similar nonuniformity of measured results. 

The flight measurements reviewed In this report seem 
generally to be quite self consistent, perhaps more consistent 
than are wind tunnel measurements. One wonders whether 
the presence of opposite walls in the wind tunnel may 
Introduce subtle local perturbations of the boundary 
layer that are only detectable as changes in the structure 
of the pressure fluctuations. Such subtle perturbations 
may strongly affect local measurements of pressure fluctuations 
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or the correlation properties of the pressure field without 
being representative of the overall boundary layer. 

Wind tunnels have been used for studying panel response. 
It is quite likely that if local perturbations to the 
boundary layer do exist as suggested above, that the overall 
response of test panels will not be seriously affected, 
since the panel may be thought of as integrating the forces 
that it sees. It thus may act to smooth or average the 
pressure fluctuation field in the same way as flight 
hardware might do. If this is Indeed the case, it would 
appear that measurements of the detailed space-time properties 
of the pressure field could not provide information that 
can readily be related to structural response. This line 
of argument leads to the recommendation that wind tunnel 
Investigations should emphasize response measurements and 
deemphasize pressure measurements. 

However, most of the response experiments performed 
to date have been concerned with displacement response 
of relatively small panels. If the linear dimensions of 
these experiments are scaled to sizes of Interest to the 
designers of high performance aircraft, the scaled frequency 
ranges fall well below the frequency ranges of interest. 

For the purpose of the designer, the models investigated 
should simulate real structures and should cover the fre- 
quency range of interest. If such programs are undertaken, 
then particular care should be taken to properly scale such 
pertinent parameters as the hydrodynamic coincidence fre- 
quency and the boundary layer thickness. Such measurements 
should be performed in conjunction with suitable flight 
measurements in order to verify the scaling relationships 
utilized In designing the model experiment. 

The data on flight measurements (ref. 39) and wind 
tunnel measurements of structure response (ref. 27) indicate 
that the use of one-to-one correspondence between aerodynamic 
and acoustic excitation can be very misleading. Also the 
use of Inflight pressure data can be misleading if the 
pressure transducer Is located near a region of strongly 
perturbed flow confined to a very local region about the 
transducer. To avoid these problems, it is recommended 
that "equivalent acoustic fields" such as presented In 
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this report be used for estimating the response of flight 
hardware and designing ground acoustic tests to simulate 
the flight environment. 
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APPENDIX A 


INVESTIGATION OF BOUNDARY LAYER PARAMETERS AT 
HIGH REYNOLDS NUMBERS AND MACH NUMBERS 

At subsonic Mach numbers and low Reynolds numbers , 
certain well defined relations exist between the boundary 
layer parameters. These relations are essentially inde- 
pendent of the Mach number, and are defined by the charac- 
teristic length and the Reynolds number. For example, 
Schlichting (ref. Al) gives the following relations for 
the boundary layer thickness 


6 

x 


0.37 Re" 1/5 » 


(Al) 


the displacement thickness 


6 


* 



(A2) 


and the momentum thickness 


0 



(A3) 


These relations follow when the boundary layer profile 
follows the well known seventh power law. At high Reynolds 
numbers and Mach numbers, the experimental evidence reviewed 
here indicates that the relationships given above must be 
considerably modified. 


Boundary Layer Thicknesses 

In figures Al, A2, and A3 we have plotted data from 
- several sources (refs. A2-A7) to determine the behavior of the 
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boundary layer thickness, the boundary layer shape parameter, 
and the ratio of boundary layer thickness to momentum thick- 
ness. The form of the presentations is the result of empiri- 
cal experimentation with the data to determine Mach number 
and Reynolds number dependence. In each figure we have indi- 
cated our choice of an empirical equation which seems both 
to fit the data best and to take on a limiting form at low 
Mach numbers and Reynolds numbers consistent with equations 
Al, A2, and A3. 

Some remarks about the data shown in the figures are in 
order. In many cases a good deal of additional work was 
required to obtain values of such parameters as displacement 
thickness from the original published data. In a few cases 
judicious guesses had to be made to fill in partially com- 
plete information. In drawing the empirical curves of best 
fit, a knowledge of the reliability of the data shown was 
kept in mind. In any case the final results, the empirical 
equations, are based on the best information available to 
us and have been derived using sound engineering principles. 
The empirical system of equations given in the figures Al, 

A2, and A3 were used to determine the following set of 
equations. 
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0.37 Re 


-1/5 
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1 + 


f „ A 
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Re , ' 
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6.9 x 10^ 



1/10 


(A4) 


6 n.3 + 0.43M g ][5/x] 

x o 8 1/3 

10.4 + 0 , 5M^ [1 + 2 X 10 Re x ] 


(A5) 


£ T75 (A6) 

X o-8 1/3 

10.4 + 0.5M^ [1 + 2x10 Re x ] 


Figures 1, 2, and 3 of the text were prepared using these 
formulas . 
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Local Coefficient of Skin Friction 

A voluminous literature exists describing the local 
coefficient of skin friction. The measurement of the 
quantity seems to contain a number of pitfalls, so that 
the calculation of the local skin friction coefficient 
from empirically based equations appears to lie in a somewhat 
controversial field. Vie have chosen to follow the lead 
of Coles (ref. A 7) and have attempted to extend his formula- 
tion to high Reynolds numbers. 

In figure A4 we have plotted the local coefficient 
of skin friction against Reynolds number based on momentum 
thickness, where both parameters are expressed in trans- 
formed coordinates based on a reference sublayer temperature 
(ref. A7) . The references from which the data were taken 
are given in the figure. The solid line represents Coles’ 
empirical approximation to all of the data available to 
him at the time of publication of his work, September 1962. 
The dotted line represents our extension to Coles’ work 
based on the plotted points. The plotted points shown in 
figure A4 show more scatter and are generally somewhat 
below Coles’ solid line and our dotted extension. The 
scatter of our data, which is greater than that of Coles’ 
reduced data, is thought to be due to the approximations 
that we have made in representing the boundary layer. 

We have used the comparison between our plotted data points 
and Coles ' line where they overlap to determine the adequacy 
of our approximations. This comparison has also been used 
to determine the placement of our dotted extension. 


Procedure for Computing c f 

The following procedure was used to enter the data 
in figure A4. It is based on an interpretation of Coles' 
work, and it contains a number of approximations to simplify 
the necessary computation. Coles shows that all available 
data on measured local skin friction for a turbulent 
boundary layer may be plotted on a single line if the 
appropriate parameters are referenced to local conditions 
at a mean position in the boundary layer. He finds that 
the local Reynolds number based on the distance from 
the wall to the mean position should have the value 

Re s = 8500. (A7) 
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This value will be used to compute the mean position 
temperature relative to the free stream static temperature. 


The following relations will be assumed to hold for 
air where the ratio of specific heats 7 = 7/5. 



1 + 0.2Mfr 
1 + 0.2M 2 r 


r = 0.89 
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, a = 0.905 


(All) 
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(A12 ) 


Equation A8 follows from equation 15.7 of reference A8 
and holds for air. It is assumed that the recovery factor r 
is in general different from unity and equal to 0.89. 

Equations A9 and A10 follow from the assumption of a perfect 
gas. Equation All holds fairly well for air over the 
temperature range from 70°K to 350°^, according to reference 
A9. Equation A12 seems to fit velocity profile data available 
to us well enough for our purpose (ref. A6*Alo). 
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We proceed as follows "by forming the relation 


Re y p P “to y 

Re 6 U„ P» u 6 


(A13) 


Use of equations A9 through A12 enables the transformation 
of equation A13 to read 
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We introduce our expression for the boundary layer thickness* 
equation A4 and write 


Re 6 - »e x | - 0.37 Re/ 5 
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(A15) 


From equation A8 we obtain 
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Use of equation Al6 allows us to write equation Al4 as 

4 


Rey - Reg 
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We set the boundary layer Reynolds number given above equal 
to Coles' value given by equation A7. The corresponding 
temperature becomes the mean sublayer temperature T . Thus, 
equation A17 becomes s 


12.31 = Re' 
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(A18) 


Equation Al8 defines the quantity Tg/T^ in terms of the free 

stream Mach number and the Reynolds number based on charac- 
teristic length. This approximate equation was used to 
compute the value of the temperature ratio which was then 
used to compute the transformed local skin friction coef- 
ficient and Reynolds number. For the latter purpose the 
following equations were used 
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Procedure for Computing c f 

The following procedure was used to calculate the curves 
for local skin friction coefficient shown in figure 4 of the 
text. The method makes use of the solid curve with its dotted 
extension shown in the figure A4 and is based on an 
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interpretation of Coles’ work due to Roshko (ref. All). The 
steps involved in the computation are as follows: 

1. Assume values of Reynolds number Re and Mach number M, 

A 

Determine 0/x using figure 3 of the text and compute Reg. 

2. Assume a value of k and compute Rbg- = k Reg, 

(One may start with a value of unity for k. ) 

3. Determine a value of c^ from Coles' curve in figure A4, 

4. Compute 


T s 

= 1 + 17.2 

( 0.022Mf ^ 

if? 

(_ 

°fc 8 

CM 

• 

o 

7* 

c f 

T w 

\1 + 0.178Mf j 

V 2 5 ^ 

l 1 

+ 0. 178mJ j 

2 


(The above equation follows from Coles’ work.) 

5. Compute 




'V ^W 1 

V T w T S 1 + 0.2M^y 


A a 


6. Compare k of step 5 with k of step 1. 

a) if they are the same, proceed to step 7 

b) if they are different, use the k of step 5 or a new 
choice of k and repeat steps 2 through 5. 


7. Compute the ratio 
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8. Compute 
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APPENDIX B 


PROCEDURE FOR ESTIMATING BOUNDARY LAYER PRESSURE SPECTRA 


The following simplified procedure may be used for 
estimating boundary layer pressure spectra for Mach numbers 
up to 4 and Reynolds numbers based on characteristic length 

from 10^ to 10 . It will be found to be of particular 
use for new vehicle configurations, when a minimum amount 
of aerodynamic Information may be available. 

Given: 

Free stream Mach number 

q Free stream dynamic pressure, in lbs/sq ft 

U^ Free stream flow velocity, in ft/sec 

x Distance from inception of turbulence (assume 

approximately 10 ft or distance from leading 
edge, whichever is greater) 

l) Calculate the overall fluctuating pressure level 

FPL overall " 20 lc «10 ■» + 84 dB < B1 > 

The standard reference sound pressure of 
2 x 10" 4 dyn/sq cm has been used.* 


* Equation (Bl) is based on the Integration of the nor- 
malized spectrum of Fig. 20 of the text. The value 
of the integral, in conjunction with a mean value 
for the local coefficient of skin friction In the 
range up to Mach 2, then gives the ratio of overall 
pressure fluctuation level to free stream dynamic 
pressure as 0.0054. This value has been used in the 
equation. It gives a value for overall sound pressure 
level which may be three decibels too high for Mach 
numbers greater than two. 
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2) Calculate the displacement thickness 


# 0.0016 x , 0 < M„ < 2 

B - (B2) 

0.004 x , 2 < < 4 


3) Calculate the characteristic frequency 

U w 

- 0.1 pp (B3) 

4) Utilize the value of f and Pig. B1 to obtain 

the estimated boundary layer pressure spectrum. 

The fluctuating pressure level in 1-Hz bands 
is obtained by adding the quantity (FPL overall * 

10 log f ) to the vertical scale of Fig. Bl, 

The spectrum shape of Pig. Bl pertains to unperturbed 
flows. Strong perturbations tend to increase 
the portion of the spectrum below f 0 » These 

increases are typically 10-20 dB at 10*" 1 f Q , 

The engineering procedure presented here is valid for 
a range of Reynolds numbers and Mach numbers which essen- 
tially cover the range of interest for present-day booster 
systems and supersonic aircraft. The effects of changing 
the aerodynamic parameters may be evaluated from the 
figures in the main body of this report. 
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FIGURE 6 





















10 log Uco F(f) - 10 log 2nq 2 c f 2 8 * dB 


+30 


CONTINUOUS TUNNEL 
60" x 84" TEST SECTION 



FIGURE 10. TUNNEL MEASUREMENTS OF BOUNDARY LAYER 

PRESSURE FLUCTUATIONS 












FIGURE 13. TUNNEL MEASUREMENTS OF BOUNDARY LAYER 

PRESSURE FLUCTUATIONS 









FIGURE 16 



TUNNEL MEASUREMENTS OF BOUNDARY LAYER 
PRESSURE FLUCTUATIONS 
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BLOWDOWN TUNNEL 
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TUNNEL MEASUREMENTS OF BOUNDARY LAYER 
PRESSURE FLUCTUATIONS 
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FIGURE 21. LONGITUDINAL SPACE CORRELATION OF THE BOUNDARY LAYER 

PRESSURE FLUCTUATIONS 
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TRANSVERSE SPACE CORRELATION OF THE BOUNDARY LAYER 
PRESSURE FLUCTUATIONS 


FIGURE 22. 
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REF. 11 
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FIGURE 23. VARIATION OF BROAD BAND CONVECTION VELOCITY 

WITH SPATIAL SEPARATION 
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FIGURE 24. TURBULENT BOUNDARY LAYER CONVECTION VELOCITY DERIVED 
FROM NARROW BAND LONGITUDINAL SPACE-TIME CORRELATION MEASUREMENTS 
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1/3 OCTAVE BAND CENTER FREQUENCIES - CYCLES PER SECOND 

FIGURE 25. TYPICAL RADIATED POWER SPECTRUM OF A THIN METAL PANEL 
EXCITED BY BOUNDARY LAYER PRESSURE FLUCTUATIONS 
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FIGURE 27. FLIGHT MEASUREMENTS OF BOUNDARY LAYER PRESSURE FLUCTUATIONS 
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FIGURE 29. FLIGHT MEASUREMENTS OF BOUNDARY LAYER PRESSURE FLUCTUATIONS 
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F I G URE 30. FLI GHT MEASUREMENTS OF BOUNDARY LAYER PRESSURE FLUCTUATIONS 
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FIGURE 31. FLIGHT MEASUREMENTS OF BOUNDARY LAYER PRESSURE FLUCTUATIONS 
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AIRPLANE NOSE CONE 



FIGURE 33. 
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FIGURE 35. FLIGHT MEASUREMENTS OF BOUNDARY LAYER PRESSURE FLUCTUATIONS 
AND AN EQUIVALENT ACOUSTIC FIELD BASED ON MEASURED STRUCTURAL RESPONSE 

(see text) 
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FIGURE 37 FLIGHT MEASUREMENTS OF BOUNDARY LAYER PRESSURE FLUCTUATIONS 
AND AN EQUIVALENT ACOUSTIC FIELD BASED ON MEASURED STRUCTURAL RESPONSE 

(see test) 
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FIGURE 39. FLIGHT MEASUREMENTS OF BOUNDARY LAYER PRESSURE FLUCTUATIONS 
AND AN EQUIVALENT ACOUSTIC FIELD BASED ON MEASURED STRUCTURAL RESPONSE 

(see text) 
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FIGURE 40. FLIGHT MEASUREMENTS OF BOUNDARY LAYER PRESSURE FLUCTUATIONS 
AND AN EQUIVALENT ACOUSTIC FIELD BASED ON MEASURED STRUCTURAL RESPONSE 

(see text) 
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